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Allgemeinste Formulierung der Ausgleichungsrechnung 
nach der Methode der kleinsten Quadratsummen 

Von Rudolf Bock, Genf (früher DGFI München) 

I. 

Die Ausgleichungsverfahren nach der Methode der kleinsten Fehlerquadratsummen sind im 
Begriffe, im Zeichen der Automatisierung von Rechenvorgängen immer größere Bedeutung 
anzunehmen. 

Während bis vor kurzem in erster Linie Geodäten und Astronomen die damit verbundene 
umfangreiche Rechenarbeit auf sich nahmen und auch einen erheblichen Anteil der Lite- 
ratur hervorbrachten, beginnen nun auch andere Teilgebiete der angewandten Mathematik 
die Ausgleichungsrechnung zur Lösung ihrer Aufgaben heranzuziehen. Die bisherigen Dar- 
stellungen der Probleme sind meist, den Gegebenheiten entsprechend,, auf die Rechnung 
mit Handrechenmaschinen abgestellt und lassen die Anwendung in modernen Rechen- 
anlagen mühsamer erscheinen, als sie in Wirklichkeit ist. 

Die Schematik der Behandlung wird vom Standpunkt der Programmierung in anderer 
Hinsicht gefordert als von dem des Formularrechners : Ein Ausgleichungsprogramm soll in 
der Lage sein, möglichst viele verschiedene Probleme durch den gleichen Satz von Routinen 
zu lösen. 

In dieser Richtung will die vorliegende Arbeit ein Schritt sein. Die Grundlagen der Me- 
thode der kleinsten Quadratsummen werden nicht erörtert (s. dazu [3], [4], [6]), ebenso- 
wenig wie die angegebenen Formeln grundlegend Neues enthalten. 

IL 

In einer Programmbibliothek darf man gewöhnlich Subroutinen für die geläufigen Ma- 
trizenoperationen erwarten. Die Ausgleichungsrechnung ist zudem für eine kurze und über- 
sichtliche Darstellung in Matrisnotation bestens geeignet. Für fast alle vorkommenden Be- 
ziehungen wird daher diese Schreibweise benutzt. Die in der europäischen Literatur übliche 
Darstellung von Matrizen durch deutsche Buchstaben (z. B. [3]) wird zugunsten der neueren 
Gewohnheit fallengelassen, Matrizen wie beliebige andere mathematische oder physikalische 
Begriffe nach Bedarf zu benennen ; die Bezeichnung von Buchstaben in halbfetter Antiqua 
ist [1] entnommen und lehnt sich an die vor kurzem in dieser Zeitschrift veröffentlichte 
Arbeit von Linkwitz zum gleichen Thema [10] an, wenn auch in der Wahl der Symbole 
letztere Veröffentlichung nicht mehr berücksichtigt werden konnte. Die Bezeichnung von 
Matrizen vom Typ (n; 1), also von einspaltigen Matrizen, als Vektoren unter Zuordnung 
kleiner Buchstaben wird beibehalten. Transponierte Matrizen werden durch Superskript T 
gekennzeichnet. "Über die benötigten einfachsten Regeln der Matrizenschreibweise infor- 
miert [3]. Näheres in [1] oder [9]. 
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III. 

Fehlerfortpflanzung. i 
Wir nehmen eine Beobachtungsreihe &i ...b n an, zu der die Gewichte <7i = ... 

j (TiCTi 

■ •• 9n = gehören (a = mittlerer oder Standard-Fehler). 

Zur linearen Funktion der bt 

f=Z a > b > (!') 

V 

gehört der Standardfehler oy, wobei 

a f a f— X a " ,a " a '' a ''- (* a ') 

V 

Nehmen wir mehrere lineare Funktionen an 

/* = 2>v&, (2') 

V 

mit den sich aus (1 a') ergebenden Standardfehlern und definieren nun eine lineare Funktion 
dieser /«zu 

9 = Y I b efe = 2>gI>!>A (3') 

Q Q » 

so erhalten wir den hierzu gehörigen Fehler nur, wenn wir die Verbindung mit der ur- 
sprünglichen Beobachtungsreihe herstellen, wie es in (3') bereits geschehen ist. Da man 
nach Umstellung erhält 

9= X b 'X b e a "Q 
v e 

ergibt sich 

a 9 a g = X a 'OAX b e a 'Q) (Z b e a ve)- ( 3a ') 

v g g 

Wir schreiben die bisher beschriebenen Vorgänge in Matrixnotation an; aus den Ge- 
wichten zum Beobachtungsvektor b wird eine (diagonale) Gewichtsmatrix G& gebildet, der 
Grund dafür ergibt sich später. Die Inverse zur Gewichtsmatrix, die nach obiger Definition 
der Gewichte in der Diagonale die Elemente OiOi stehen hat, nennen wir Korrelations- oder 
Fehlermatrix 6?6 -1 . Im speziellen Fall der 6, wo die Korrelationsmatrix wieder Diagonal- 
form annimmt, sprechen wir von unkorrelierten oder statistisch unabhängigen Beobach- 
tungen. Wir verwenden dieselbe Gleichungsnumerierung wie oben und erhalten aus 

/=arb (1) 

(la) 

(2) 

(2a) 

wobei der Typ der Matrizen nun allgemein gehalten ist. 
Weiter gilt für eine Funktion von i aus (2) 

g = Bf=BAb (3) 

die Korrelationsmatrix 

Gg 1 = (BA) Gjj^BAjT = BAG^A^B^ = BGj 1 «*'. (3a) 

womit nun die Bedeutung der quadratischen Korrelationsmatrix erbracht ist: (2a) sagt 
über die Fehler der Komponenten von f mehr aus als (2a'), welches nur die Diagonal- 
elemente von Gj 1 in (2a) beschreibt. Die Kenntnis der gesamten Korrelationsmatrix ist 
gerade ausreichend, um die Fehler linearer Funktionen der f zu bestimmen, ohne auf die 
ursprünglichen, nichtkorrelierten Beobachtungen zurückzugehen. 



für das Fehlerquadrat 






ö / -i = a^G6 1 a. 


Die Erweiterung auf 






i = Ab 


liefert 






Gj x = AG; 1 Ar, 
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Wir fassen zusammen : 

Die lineare Transformation eines Vektors 

x = My (4) 

zieht eine entsprechende Umformung der zugehörigen Fehlermatrix nach sich : 

Gi 1 = MG; 1 Mr. (4a) 

Es ist noch wichtig, auf die Bedeutung der Elemente der G _1 hinzuweisen. Aus der 
Statistik (s. [2]) ergibt sich, daß das Element ik einer solchen Matrix den mittleren Wert 
des Produktes ff* ff* (= Kovarianz) angibt, wenn Oi und o* die dem i-ten und fc-ten Element 
des untersuchten Vektors zugeordneten mittleren Fehler sind. Daß ein Diagonalelement 
gleich dem Quadrat des entsprechenden mittleren Fehlers („Varianz") ist, ergibt sich daraus 
von selbst. Beachten wir aber, daß wir, um diese Bedeutung zu wahren, an keiner Stelle 
einen Skalarfaktor zu Gewichts- oder Fehlermatrizen zulassen dürfen, sondern die zu 
Anfang eingeführte Definition von Gewichten und Elementen der Fehlermatrizen bei- 
behalten müssen. Siehe hierzu auch [10], S. 160 ff. 

Man schreibt häufig, um den Mittelwert des Produktes zweier Fehler anzudeuten, die 
Form CTä~ö*- 

Die Gewichtsmatrizen spielen eine nicht so leicht überschaubare Rolle wie die Fehler- 
matrizen. Wir definieren sie lediglich als Inverse der letzteren und werden später sehen, daß 
damit tatsächlich teilweise eine formale Verwendung als Gewichtsmatrix möglich ist. Vor- 
aussetzung ist natürlich, daß G _1 nichtsingulär ist. Siehe hierzu die Anmerkungen am Ende 
von Abschnitt VII. 

IV. 

Im Fall einer nichtlinearen Funktion ist das Vorgehen analog (4) und (4a), wenn nur 
vorher durch Taylorentwicklung eine linearisierte Form an der entsprechenden Stelle her- 
gestellt wurde. 

Um auch die Differentiation unter bestimmten Voraussetzungen in Matrizenschreibweise 
durchführen zu können, definieren wir die Ableitung einer Funktion eines Vektors nach 
diesem Vektor als diejenige Matrix, die als Element ik die Ableitung der Komponente i der 
Funktion nach der Komponente k des Vektors enthält. Es ist leicht zu zeigen, daß im Fall 
einer linearen Funktion die Ableitung die Form annimmt 

dx dx 

Nur in dieser Definition lassen wir Differentiationen nach Vektoren zu. 

Für den Fall der Fehlerfortpflanzung bei nichtlinearer Funktion erhalten wir nun, wenn 



die linearisierte Darstellung 



und daraus 



f = f(x) (6) 

dl 

f(x + dx) = f(x) + — dx (6a) 

dx 

*■-(£)<<)'• 

Wir sehen, daß (4a) lediglich einen Sonderfall von (6b) für eine lineare Abhängigkeit f(x) 
darstellt. 

(wird fortgesetzt) 
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Allgemeinste Formulierung der Ausglekhungsrechnung 
nach der Methode der kleinsten Quadratsummen 

Von Rudolf Bock, Genf (früher DGFI München) 
(Forlsetzung von Seite 46) 



Bevor wir den allgemeinsten Fall in Angriff nehmen, behandeln wir kurz die beiden 
klassischen Fälle der Ausgleichungsrechnung und beginnen mit dem Problem der ^Aus- 
gleichung nach vermittelnden Beobachtungen". 

Hierbei liegt bekanntlich ein Beobachtungsvektor b mit bekannter (diagonaler) Korrela- 
tionsmatrix Gjj 1 vor und soll einer (zunächst linearen) Funktion eines Vektors x von Un- 
bekannten gleichgesetzt werden. Diese Funktion drücken wir aus durch eine Matrix A und 
ein konstantes Additionsglied c. 

Die Gleichsetzung kann bei Überbestimmung nur erfolgen, wenn die Beobachtungen 
durch Verbesserungen v korrigiert werden und die „Fehlergleichungen" lauten: 

Ax + c + b + v = 0. (7) 

Dabei sollen die x so bestimmt werden, daß 

M = v T GbV = Minimum. (7a) 

Aus (7) nehmen wir v = — (Ax + b + c) und setzen iii (7a) ein: 

M = (Ax + b+ c)^Gb(Ax + b + c). (8) 

Die Erfüllung von (7a) führen wir herbei, wenn wir M einzeln nach den Komponenten 
von x partiell differenzieren und jeden dieser Differentialquotienten gleich Null setzen. Nach 
unserer in TV. eingeführten Definition können wir schreiben 

£-.. 

Die eingeführte Ableitungsregel ist anwendbar, obwohl in M ein x T auftritt : Wir differen- 
zieren zunächst nach x und betrachten x T als zweite Variable. Nach dieser differenzieren 
wir, wenn wir M durch M T ersetzt haben, welche zwei Größen natürlich identisch sind. 
Diese zweite Differentiation ist dann der ersten gleichlautend und das Ergebnis ist (nach 
Transposition, um eine Vektorgleichung zu erhalten): 

(^) T =2ArG b (Ax + b + c) = (!•) 
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womit die bekannten Normalgleichungen beschrieben sind. Schematisch lautet die Lösung 

x= -(ArG b A)-iArG„(b+c) (11) 

was, in (7) und (8) eingesetzt, auch die v und M zu berechnen gestattet. 
Auch die als Rechenproben viel benutzten Beziehungen 

M = v^GbV = (b + c)^G b (Ax + b + c) = - (b + c)^G b v (12) 

leiten sich aus (7a), (8) und (10) ohne weiteres ab. Sie entsprechen in Summenschreibweise 
den Formen 

M = [gvv] = [gll] - [Xilgaj]] = - [glv] (12') 

wobei [gll] — [arjfsra^]] = [gll ■ n], denn die Produktsumme aus Lösungen und rechten 
Seiten der Normalgleichungen ist gleich dem Reduktionsbetrag, der bei der üblichen Gauß- 
schen Elimination an [gll] angebracht wird. 

Es bleibt noch übrig, die sich aus (11) abzulesende lineare Abhängigkeit zwischen den x 
und b zu benutzen, um eine Korrelationsmatrix für die x abzuleiten : 

G; 1 = [(A^G fc A)-iArG b ]Gi; 1 [(ArGbA)-iArG b ]y = 

= (ArGbAJ-iArGbG^GbA^GbA)-! = (13) 

= (A^G b A)-i. 

Auch der Funktion b -)- v ist eine Korrelationsmatrix zugeordnet. 

Wir stellen b + v als lineare Funktion der b dar, was aus (7) und (11) ergibt: 

b + v = - (Ax + c) = A(A r G b A)-iA^Gb(b + c) - c (14a) 

G( b 1 +¥) = A(A^G b A)-iA^. (14) 

Die Quadrate der mittleren Fehler der ausgeglichenen Beobachtungen bzw. der ermit- 
telten Unbekannten erhält man in bekannter Weise aus den Diagonalelementen der Kor- 
relationsmatrizen durch Multiplikation mit dem Faktor M /(»ib — «i) (wb = Zahl der Be- 
obachtungen, n z = Zahl der Unbekannten). Sind die eingeführten mittleren Fehler der I» 
statistisch richtig gewählt, so ist ein mittleres M = «b — n* zu erwarten (% 2 - Verteilung). 
Der Rang der beiden Matrizen Gj 1 und G^J" +V ) ist offensichtlich höchstens der von A, näm- 
lich n x . Siehe hierzu die Anmerkungen im Anschluß an (35). 

Für den Fall nichtlinearer Abhängigkeit 

l'(x) + b-fv = o (7b) 

erhält man einen Ansatz (7) mit c = f (x), A = df/dx wenn x eine Näherung für x ist, die 
auf beliebigem Weg vorherbestimmt ist. 

An Stelle von x bestimmt man nun ein Ax v und es ist x»+i = x, -f- Zlx, + i. Der Index gibt 
dabei die Iterationszahl an, da die Bestimmung von Zlx mehrmals wiederholt werden muß, 
um der Änderung von A mit x Rechnung zu tragen. Über Ausgangskriteria des Iteration.s- 
zyklus siehe Abschnitt VII. 

VI. 

Wir gehen über zum zweiten Standardfall der Geodäsie, der „Ausgleichung nach be- 
dingten Beobachtungen". 

Die Bedingungsgleichungen seien formuliert durch 

B(b + v) + c = (15) 

wobei B die spezielle Form der Bedingungen charakterisiert und die übrigen Buchstaben 
die Bedeutung wie in V. haben. G^ 1 wird wieder als bekannt vorausgesetzt. Die Minimums- 
bedingung ist die gleiche wie oben: 

M — v T GbV = Minimum . (7a) 
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Bekanntlich löst man das Problem durch Einführen der Lagrangeschen Multiplikatoren a 
und erweitert (9) in folgender Form 

M = vrG„v + 2ar(B(b + v) + c). (16) 

Die partiellen Ableitungen von M nach den Elementen von v müssen Null gesetzt werden 
und nach (5) erhalten wir: 

(^f=2G b V + 2BTa = 0. (17) 

Daraus folgt 

v=-Gi 1 B^a (17a) 

und über (15) 

a = (BGi 1 Br)-i(Bb + c), (18) 

wobei die Elemente von Bb + C gewöhnlich als Widersprüche bezeichnet werden. 

Die Lösung, nämlich v, ist bereits durch (18) und (17 a) bestimmt. Wir schreiben die 
Korrelationsmatrix der a an, die sich nach (18) ergibt: 

G; 1 = [(BGü 1 Br)-iB]Gi 1 [(BGi; 1 Br)-iB]2'= (BG^B^-i. (19) 

Aus (17a) ergibt sich weiter 

M = v^GbV = a^BG^Bra = a^G.a (19a) 

und als Rechenprobe noch 

M = \Tq h y = a r(Bb-r c) 

(19 a) zeigt, daß für den Fall der a die Einführung der Inversen zur Korrelationsmatrix als 
Gewichtsmatrix gerechtfertigt ist (siehe hierzu [10]). 

Aus einer Darstellung von b + v als lineare Funktoin von b 

b + v = b-G^B^a^ (E-Gb 1 Br(BGi 1 Br)-iB)b + const. (20a) 

(E = Einheitsmatrix) 

ergibt sich die Korrelationsmatrix der b + v 

G£ +y) = Gü 1 - G^B^BG-jB^-iBGi 1 . (20) 

Aus BG^ +T j = ergibt sich, daß G("^ +T ) singulär ist mit einem Rangdefekt, der gleich ist 
der Anzahl der Zeilen von B, also gleich der Anzahl ni der Bedingungsgleichungen (15). 
Siehe hierzu die Anmerkungen zu Gleichung (35). 

Die Anzahl der Freiheitsgrade für M ist ebenfalls gleich n\. 

Im Fall nichtlinearer Bedingungsgleichungen f (b + v) + c = stellt man die Linear- 
form durch Differentiation an der Stelle b her und erhält 

rlf 

Bv + f(b) + c=0 mit B = — . (21) 

db 

Sind die für eine passende Lösung gestellten Kriteria (s. Abschnitt VII.) jedoch streng, 
so genügt ein einmaliges Durchrechnen mit diesem Ansatz nicht, da ja die v nicht wirklich 
differentiell klein sind. Auch die Lösung des Problems der bedingten Beobachtungen er- 
fordert bei nichtlinearen Bedingungsgleichungen in Strenge ein iteratives Vorgehen. Man 
schreibt 

Bv„ +1 + f(b + v,) + c-Bv,, = (22) 

mit 

df. 
dbJk=k +Tr 

Damit ist eine (15) entsprechende Form gefunden. An die Stelle von Bb + c tritt 
f (b + v») + C — Bv,,. Diese Größe, die dem ursprünglichen Widerspruch gleichsteht, geht 
mit wachsender Iterationszahl nicht gegen Null, sondern gegen den Wert — Bv. 



B = 
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VII. 

Wir formulieren nun die allgemeinste Form des Ausgleichungsproblems. Ausgehend von 
einem Satz von Bedingungsgleichungen zwischen Unbekannten dx, Beobachtungen b und 
Verbesserungen v schreiben wir direkt die Linearform an : 

f = Az1x + Bv+c = (23) 

wobei wir bereits wissen, daß c = c(b) und daß im Fall nichtlinearer Bedingungsgleichun- 
gen f (x, b + v) = die Matrizen die Bedeutung haben 



A = 



df 



dx 



und B = 



df 



db 



b = b + T, 



und außerdem c = c(x„, b + v„) (vgl. (7b) und (21)). 

Zusätzlich zu (23) ist wieder M = v^GbV zum Minimum zu machen. Die Ähnlichkeit 
der Aufgabe mit der Formulierung (7) ist auffallend, lediglich ist v durch Bv ersetzt. Daß 
rein formal nur G^ 1 und Gb durch die aus (4) und (4a) folgenden Matrizen Gb 1 = B G^B 7 
und Gb = (BG^B 7 ) -1 zu ersetzen sind, bedarf jedoch erst der Ableitung. Wir führen die 
Lagrangeschen Multiplikationen ein und schreiben 

M = vrG b v + 2ar(AZlx + Bv + c). (24) 

Die Differentiation nach v liefert (s. auch (17)), 

AM\T 

— =2Gi,v + 2BTa=0 (25) 

und daraus wieder (17 a) 

v=_G b 1 Bra. (26) 

Dieses Ergebnis setzen wir in (23) ein und erhalten 

a = (BGi 1 Br)-i(Azlx + c) = G„(AZlx+ c) (27) 

Die Differentiation von M nach Zlx liefert ausgehend von (24) 

AM 

T - r = 2a^A=0 (28) 

dzlx 

liefert zusammen mit (27) die Bestimmungsgleichung für Zlx: 

2A T Gu(AZlx + c) = (30) 

woraus sich die Lösung 

Ax = - (A^GBAJ-iA^GBC (31) 

tatsächlich in der oben erwähnten Analogie zum Fall der „vermittelnden Beobachtungen" 
ergibt. Auch für M läßt sich schreiben 

M = (Bv)rG B (Bv). 

Als Anzahl der Freiheitsgrade erhalten wir %i — n z . 

Die unter V. und VI. beschriebenen Ausgleichungsverfahren lassen sich nun ohne Schwie- 
rigkeiten in das eben erhaltene Schema einordnen. Sämtliche dort erhaltenen Ergebnisse 
erhält man auch aus (23) bis (31), wenn für die Matrizen die Spezialformen eingeführt 
werden: Für den Fall der vermittelnden Beobachtungen B = E (Einheitsmatrix), für den 
Fall der bedingten Beobachtungen ein A vom Typ (n;0), d.h. es existieren keine Un- 
bekannten. Auch die Zahl der Freiheitsgrade resultiert zwanglos aus diesen Vereinfachungen 
des allgemeinen Falles: Für B = E ist die Zahl der Gleichungen gleich der (gewöhnlich 
irrelevanten) Zahl der Beobachtungen, im anderen Fall ist die Anzahl der Unbekannten 
Null. Aus den Fehlergleichungen wird bei dieser Betrachtungsweise lediglich eine spezielle 
Form von Bedingungsgleichungen. 
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Wir betrachten noch kurz die Lösbarkeitsbedingungen des Ansatzes (23). Die Normal- 
gleichungen für <x enthalten Skalarprodukte von Zeilen aus B und BGj; 1 . Die Zeilen müssen 
also linear unabhängig sein und dürfen nicht gleichzeitig lauter Nullelemente enthalten. 
Überdies muß B wenigstens dieselbe Anzahl von Spalten wie Zeilen haben, d. h. die Zahl 
der voneinander unabhängigen Beobachtungen muß wenigstens gleich sein der Zahl von 
Bedingungsgleichungen*). Gleiches wie für die Zeilen von B ergibt sich aus den Normal- 
gleichungen der Ax für die Spalten von A. Die eingeführten Unbekannten müssen un- 
abhängig und ihre Zahl darf höchstens gleich der Zahl der Bedingungsgleichungen sein. 

Für die ausgeglichenen Werte sind noch die zugehörigen Korrelationsmatrizen abzuleiten. 
Hierzu stellen wir wieder die Ax und b + v als Funktionen der b dar und wenden (4), (4a) 
an. Zunächst erhalten wir aus (31) und der Definition von c 

~ = - (A^GbAJ-iA^Gb ^ = - (A^CA)- 1 Ai-GiB (32) 

du ab 



woraus resultiert 



Weiter ergibt (26) bzw. (27) 



G2i=(A^G B A)-i. (33) 



d(b + v) =E+ dv_ E _ G . lBTGB ^_ A dJx )= (34) 



db ' db " \db dk 

= E - G^B^GbB + G„ 1 BrG B A(A 7 "GiA)-i A r G E B 

woraus nach einiger Matrizenalgebra folgt 

G <b + v>= Gü 1 - «VBtGbBG,; 1 + G i ; 1 BrGBA(A'-G,AriArG,BGi 1 . (35) 

Die abgeleiteten Korrelationsmatrizen (33) und (35) entsprechen bei Einführung der oben 
erwähnten Sonderfälle für die klassischen Fälle der Ausgleiehungsrechnung den Ausdrücken 
(13) und (14) bzw. (20). Führt man übrigens ein nichtsinguläres A quadratischen Typs ein, 
d. h. ebensoviele Unbekannte wie Bedingungsgleichungen vorhanden sind, so ergibt sich 
v = und G(jJ +y ) = Gjj 1 , wie zu erwarten ist. Es ist nicht so leicht wie in den Abschnitten V. 
und VI., zu zeigen, daß auch hier G^ + v) singulär ist mit einem Rangdefekt, der genau gleich 
ist der Anzahl der Freiheitsgrade. Um die dort gemachten kurzen Andeutungen zu einem 
allgemeinen Beweis zu erweitern, müssen wir etwas ausführlicher auf das Problem eingehen. 
Mathematische Hilfsmittel sind wieder [1] und [9]. Zunächst ersetzen wir den Rang von 
G(b+v) durch den von 

n ß -l r d(b + v) 

der dem von G^ +r) gleich ist (Gb ist nicht singulär). Aus D 2 = D ergibt sich, daß alle 
Eigenwerte Xt von D Null oder Eins sind, da für die Eigenvektoren yj die allgemeine 
Beziehung D 2 y s = Aj 2 y< gilt. Nun ist die Summe der Eigenwerte einer Matrix gleich ihrer 
Spur; in unserem Fall auch gleich dem Rang. Da D = Di — D2 + D3 (Gleichung (34)) 
können wir für jedes D* die Spur getrennt ermitteln und erhalten Sp(D) = Sp(Di) — 
— Sp(J)2) + Sp(J)3). Außerdem gilt für jedes i wieder (D<) 2 = D(. Wir gehen also zurück 
auf die Ermittlung des Ranges der D* und erhalten: 

Rang (Di) = Rang (E) = m, 

Rang (D 2 ) = Rang(Gb 1 ßrG B B) = n, 



*) Bedingungsgleichungen nur zwischen Unbekannten verlangen also eine Sonderbehandlung. Hierzu 
folgt man vorteilhaft der Methode in zwei Schritten nach [7] und spaltet diese Bedingungen der Form 
AAx + c = ab. Man berechnet nach (31) ein Ax und erweitert dann die „rechte Seite" — \. T ii^C von 
(30) um — A T [A(A t GbA) _1 A t ] -1 (Az1x 4- c). Die nochmalige Auflösung für zlx ergibt die endgültigen 
Werte. Vergleiche die ausführliche Darstellung in [10]. 
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(da wegen Rang (xy) ^ min (Rang (x), Rang (y) ) Wf der maximale Rang für D 2 ist, anderer- 
seits durch BD2 = B genau n\ Eigenwerte zu 1 bestimmt sind und wj daher auch Minimal- 
rang ist. 

Rang (D 3 ) = Rang (Gi 1 BrG B A(ArG B A)-iArG B B) = n % 

(da wie oben der Rang höchstens n x sein kann und durch A T G B BD3 = A 2 'G B B wieder 
gerade n x Eigenwerte = 1 gegeben sind). 
Unser Gesamtergebnis ist also 

Rang (G^ +v) ) = n h — n t + n T . 

Von Interesse sind in manchen Fällen auch die Korrelationen zwischen den Ax und b + v, 
also mit (33) und (35) die Korrelationsmatrix für alle ausgeglichenen Größen. Man erhält 
sie anschaulich, indem man aus den Elementen von b + v und Ax einen Vektor aufbaut 
und die Beziehung mittels (32) und (34) auf die b zurückführt. 

Man erhält allgemein für zwei Vektoren x, y 

dx \ n-i / d y 



und in unserem Fall 

K(k+Y.4i) = -Gi 1 B2'G B A(ArG B A)-i. (37) 

Die gesamte oben erwähnte Korrelationsmatrix hat das Aussehen 






\ A (b+T,JX) "Ji / 

Die Matrizen (33), (35) und (37) scheinen auf den ersten Blick äußerst lange Ausdrücke 
zu sein, sind jedoch aus wenigen Teilmatrizen aufgebaut, die schon im Verlauf der Rech- 
nung auftreten. Zudem werden ja Korrelationen erst nach Abschluß und Prüfung der ge- 
samten Ausgleichung berechnet und keinesfalls innerhalb des Iterationszyklus. 

VI IT. 

Das eben gegebene Stichwort „Prüfung der Ausgleichung" bedarf noch einer Erklärung; 
während man sich bei Handrechnung gewöhnlich mit schrittweisen Rechenkontrollen und 
der M-Probe analog (12) bzw. (12') oder (19a) zufrieden gibt und bei iterativer Verbesserung 
die Änderungen von Schritt zu Schritt beurteilt, kann ein automatisch arbeitendes Pro- 
gramm jedenfalls noch tiefergreifende Prüfungen der Ergebnisse enthalten. 

Handelt es sich um lineare Bedingungsgleichungen, so sind A und B von Ax und v unab- 
hängig und die Prüfung hat sich prinzipiell tatsächlich nur auf richtige Rechnung im ersten 
Schritt zu erstrecken. Die M-Probe kann allgemein formuliert werden 

l/ = v r G b v = a r (Az1x + c) (39) 

und gibt eine gute Kontrolle für grobe Rechenfehler und auch für schlecht konditionierte 
Normalgleichungen. Es empfiehlt sich, diesen von der Linearitätsbedingung unabhängigen 
Test auch bei iterativer Anwendung der Formeln durchzuführen. Ein gutes Kriterium ist 
ein Zusammenstimmen der beiden Werte auf 10 -4 (relativ), wobei jedoch für den Fall M = 
(quadratisches nichtsinguläres A) eigens ein Absolutwert-Test vorgeschaltet werden muß. 

Im Fall der Nichtlinearität verprobt man auf diese einfache Weise zwar den einzelnen 
Iterationsschritt, es bleibt aber die Frage, ob die eben abgeschlossene Iteration ein be- 
friedigendes Ergebnis gebracht hat. Wir suchen also nach Ausgangskriterien für den Ttcra- 
tionszyklus. 

Eines der wichtigsten ist zweifellos die Erfüllung der Bedingungsgleichungen, die ja je- 
weils nur in linearisierter Form zur Erfüllung gebracht werden. Der Betrag des Residuen- 
vektors |f | (= Betrag der noch vorhandenen Widersprüche) ist nur in besonderen Fällen 
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eine sinnvolle Testgröße, da die Bedingungsgleichungen ja ganz verschiedenen Ursprungs 
sein können und dann zeilenweise völlig andere Gewichte zugeordnet werden müssen. Be- 
stimmen wir eine Gewichtsmatrix aus 



*l\ «- mi 



fr 1 

db/ Ub Vdb/ 



-1 



= (BG l | 1 Br)-i = GB (40) 



so gibt <5f= {f T G B f} 1/2 (40a) 

eine erste brauchbare Testgröße, für die bei richtiger Fehlerzuordnung eine Grenze 10~ 4 
eingeführt werden kann. Ein sinnvoller Ansatz vorausgesetzt konvergiert <5f sehr rasch 
und läßt sich bis zur Rechengenauigkeit herunterdrücken. 
Mehr statistische Bedeutung hat die Korrelationsmatrix 

/ df_ dl\ / dl d M r =(B4) (C,K\H (41) 

(angeschrieben als Matrizen von Matrizen). Dabei sind G^J +T ), K und G^ x nach (33), (35) 
und (37) von der zuletzt durchgeführten Iteration v zu nehmen, B und A in (41) müssen 
jedoch bereits an der Stelle b -f- v,,+i, x,+i abgeleitet sein. Andernfalls ergäbe sich Gj" 1 s 0. 
Die Diagonalglieder von G, 1 geben die <7<<7j für die Komponenten von f nach dem Schritt v 
entsprechend der Definition (4a) f. Die Untersuchung von G^ 1 nach jedem Schritt würde 
jedoch die Rechenzeit auch bei schnellsten Anlagen erheblich verlängern. — Allein mit der 
Erfüllung der Bedingungsgleichungen kann sich eine strenge Prüfung noch nicht zufrieden 
geben, vielmehr bedarf auch die zweite in die Rechnung eingeführte Bedingung (M = Mini- 
mum) der Prüfung. Wir prüfen also die Ableitungen von M nach b + v und x, die bei der 
endgültigen Lösung verschwinden müssen. 
Wir definieren 



fV/ dM - \( dM - \1 1/2 



(41a) 

dM _ \ldM_ \1 1/2 



i; 

wobei i,j Elemente aus den entsprechenden Vektoren bezeichnen und aä = aä nach (4a) f. 
Eigentlich wäre es nötig, auch noch die Korrelationen zwischen den Ax und b + v zu berück- 
sichtigen; jedoch ist bereits die exakte Durchführung von (41a) langwierig und man be- 
schränkt sich vorteilhaft sogar auf eine Näherung für 6M± und schreibt 

1/ dilf \ , / dM \r|i/2 

(41b) 
dM\ „_, (dMY) in 



«-H-dTj^Vdx-; 

wobei G^ +T) durch Gjj 1 ersetzt wurde und damit das Kriterium strenger wird, da 

Gl — Gf b+T) = G; 
und somit positiv definit ist. 

Führen wir die Differentiationen durch, so erhalten wir entsprechend (25) 

dM \ T fdM\ T 

-1 =2G b v + 2Br a (B = B, +i !) 



^d(b + v)/ \dv ) 

und mit (26) 

dM V 

= 2(B,+i-B,)r a (42 ) 



was keinen besonderen Rechenaufwand verlangt. 
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Ähnlich benutzen wir (30) und erhalten 

lirf = ((ÖxT = 2ATG *( AAx + c ) = - 2A^G B Bv. (43) 

Auch für dM\ und &M2 oder 6M = yöMi 2 -\- 6M<p lassen sich bei richtiger Anfangs- 
annahme G^ 1 numerische Werte vorschreiben, die jedoch vorteilhaft nicht kleiner als 10 -4 
liegen, da sonst zu viele unnütze Iterationen durchgeführt werden. 

Die drei Größen öl , 6Mi und ÖJÜ2 prüfen ausreichend den zuletzt durchgeführten Schritt 
auf ein befriedigendes Ergebnis. Bei Nichterfüllung einer der gesetzten Bedingungen geht 
man in den nächsten Schritt, zählt aber die Anzahl von Schritten und läßt bei einer vor- 
geschriebenen Höchstzahl abbrechen (mit den vorgeschlagenen numerischen Werten erhält 
man bei sinnvollen Ausgangswerten selten mehr als 3 — 4 Iterationen). Um jedoch diesen 
Fall der Unlösbarkeit gegebenenfalls abzukürzen und gleichzeitig weitere Rechenkontrollen 
einzubauen, kann die Konvergenz von dl und 6M geprüft werden, d. h. jeder dieser Werte 
muß von Schritt zu Schritt kleiner werden. Ein sofortiges Abbrechen bei einmaliger Nicht- 
Konvergenz empfiehlt sich jedoch nicht, da diese oft scheinbar ist und numerische Ursachen 
hat. Für die Größe M ist eine Konvergenzprüfung unsinnig, solange die Bedingungs- 
gleichungen nicht hinreichend erfüllt sind. Numerische Beispiele zeigten, daß sich selbst 
nach Erfüllen der Bedingungsgleichungen entsprechend dem gegebenen öf-Zahlenwert M 
noch nicht genügend gut entlang der Kurve bewegte, deren Minimum gesucht ist. Es emp- 
fiehlt sich, 6M -Werte und ÖM -Konvergenz erst nach Erfüllung der Bedingungsgleichungen 
zu berechnen und zu testen. 

IX. 

Abschließend sollen die nötigen Formeln zur Lösung des allgemeinsten Falles der Aus- 
gleichungsrechnung in geschlossener Form angeschrieben werden. 

Wir beginnen mit Beobachtungen b, einer Korrelationsmatrix GjJ 1 , einer ersten Nähe- 
rung v v= o = 0, Näherungen für Unbekannte x U; einer allgemeinen Funktion f und deren 
Ableitungen in Matrizenform nach x und b. Diese Ableitung kann u. U. numerisch erfolgen, 
d. h. in die dem Ausgleichungsprogramm zugehörige Routine zur Berechnung von f wird 
mit jeweils einer geänderten Komponente x^ + Axi oder b( = y v i + ^Ibj eingegangen und 
eine Spalte von A oder B aus (ljAxt) AI oder (1/Zlb<) AI erhalten. Eventuell durchläuft 
man diesen Zyklus zweimal und ändert dazwischen die Zeichen der Axt, Abt; die zwei 
Spalten werden dann gemittelt. Die Größen A , für die f als linear betrachtet werden kann, 
müssen im Fall der numerischen Differentiation als Ausgangsinformation erwartet werden. 
Die analytische Berechnung von A und B ist zweifellos schneller. Siehe [5] und [8] zur 
numerischen Differentiation. 

Für einen Iterationsschritt haben wir die Rechenschritte (ohne Rücksicht auf Wieder- 
verwendbarkeit von Zwischenergebnissen) : 

Lösungen 

f = f(i„b + v.) 

\dx/x=x„ 

\ db/b=b+v„ 

c = f-Bv, 
Gi 1 = BGj; 1 Br 

Ax, +1 = - (A^GbA)-i A^Gbc 
a = GbAzIx + Gbc 
v, +1 = -Gi 1 Bra. 
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Fehlerrechnung 



M r . 



Berichte und Mitteilungen 



VJT +1 GbV„ + l 



zfv 

Nr. 3/1961 



(A^GbA)-i 



G^ +v )= G; 1 - G^B^GbBGj; 1 + G^' n^GitAfArGBAl-iArGBBGb 1 

K ( bi-v,4x)= -Gb ) B? , (i ll A(A'? , GBA)-i. 



Prüfung 



v T G b v - ol t (A Ax + e)\ g £l 



a) Iterationsschritt 
1 

M 

SU - öh+i ^ 
(ÖM, - ÖM V+X ^ 0) 

b) Abschluß 

<5f= {i r G B *} 1/2 g£2 
ÖMi = 2 {dfG^di} 1 ' 2 mit di = (H„ +1 - B„) T a 
< 5M2 = 2{d| , (A^G B A)- 1 d 2 } 1/;! mit d 2 = A y (i B Bv 
ÖM = {d.¥i 2 + ditf., 2 } 1 ' 2 ?S £ 3 . 
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